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Abstract 

The program of quantizing the gravitational field with the help 
of affine field variables is continued. For completeness, a review of 
the selection criteria that singles out the affine fields, the alternative 
treatment of constraints, and the choice of the initial (before imposi- 
tion of the constraints) ultralocal representation of the field operators 
is initially presented. As analogous examples demonstrate, the intro- 
duction and enforcement of the gravitational constraints will cause 
sufficient changes in the operator representations so that all vestiges 
of the initial ultralocal field operator representation disappear. To 
achieve this introduction and enforcement of the constraints, a well 
characterized phase space functional integral representation for the 
reproducing kernel of a suitably regularized physical Hilbert space is 
developed and extensively analyzed. 
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I. 



INTRODUCTION AND THE 
MAIN POINTS OF PAPER I 



In a previous paper (hereafter referred to as P-I) an introduction and 
outline of a program of noncanonical quantization of the gravitational field 
was presented. The key concepts in the present approach are (i) a careful 
selection of the basic kinematical variables, (ii) the use of a quantization pro- 
cedure that treats all constraints alike, (iii) the use of ultralocal field operator 
representations prior to introducing constraints, and (iv) the imposition of 
the gravitational constraints, a process in which all traces of the temporary 
ultralocal representation characteristics are replaced with the physically rel- 
evant ones. In this section, we briefly review topics (i), (ii), and (iii) which 
have largely been discussed already in P-I. 

A. Nature of the basic gravitational variables 

One of the central requirements for the present program is the preservation, 
on quantization, of the positive-definite character of the spatial part of the 
classical metric gab{x), a,b E {1, 2, 3} (or more generally a,b E {1, . . . , s} in 
an s-dimensional space, s > 1). Insisting on this requirement leads us to 
adopt affine commutation relations (in contrast, for example, to canonical 
commutation relations) which are expressed in terms of the local metric field 
operators gab{x) [= gba{x)] (which were denoted by aab{x) in P-I) and suitable 
local "scale" field operators 7r^(a;) (which were denoted by k^(x) in P-I). The 
word "local" here is intended to mean that these expressions only become 
operators after smearing with suitable spatial test functions. In units where 
h = 1, which are commonly assumed throughout this paper, the basic set of 
affine commutation relations reads 



These commutation relations are translations of identical Poisson brackets 
(modulo ih, of course) for corresponding classical fields, namely, the spa- 
tial metric gab{x) and the mixed- valence ("scale") field vr^(x) = gbd{x)'7i'"^{x), 
along with the usual Poisson brackets between the metric field gab{x) and 




(1) 



[gabix), gcdiy)] = . 
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the canonical momentum field 7r^'^(x). While classically there is essentially 
nothing to be gained by using the field vr^(x) rather than 7r'^'^(x), quantum me- 
chanically the situation changes completely. This change arises because the 
affine commutation relations — and only the affine commutation relations — 
admit local self-adjoint operator solutions for both gab{.x) and 7r^(a;) which 
in addition have the property that gab{x) > for all x. This latter property 
means that for any nonvanishing set of real numbers and any nonvan- 
ishing, nonnegative test function, /(x) > 0, that 



Other choices of basic variables fail this test. For example, self-adjoint canon- 
ical variables lead to metrics that have spectra unbounded below as well as 
above, while triad fields and their canonical partners lead to metrics that are 
nonnegative, but not necessarily positive definite. 

B. How quantum constraints are to be imposed 

Since gravity is a reparametrization invariant theory, it follows that the 
dynamics — indeed, the entire physical content of the gravitational field — 
enters through imposition of the relevant constraints, specifically, the diffeo- 
morphism and Hamiltonian constraints 0, |^. These constraints lead to an 
open set of first-class constraints in the classical theory, but on quantization 
and in virtue of an anomaly (or, alternatively, a factor ordering problem), 
they give rise to a set of operator constraints that, partially at least, are sec- 
ond class in nature. There exist several methods to deal with the quantum 
theory of second-class constraints in the literature, but generally such meth- 
ods treat the first- and second-class constraints, and thereby the variables on 
which they depend, in fundamentally different ways. 

An exception to the rule of a different operator treatment for first- and 
second-class constraints is offered by the projection operator method approach 
to the quantization of systems with constraints [^, ^, |^. The projection oper- 
ator method was already discussed in P-I, and that discussion also included 
several elementary applications of the technique to simple, few degree-of- 
freedom systems. We do not repeat that discussion here. Instead, we simply 
observe that, rather than impose the self-adjoint quantum constraints in the 
idealized (Dirac) form ^a\4')phys = 0, a G {1,...,A}, on vectors \4')phys 




(2) 
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in a putative physical Hilbert space, \4')phys £ ^phys, we define a (possibly 
regularized) Sjphys = 1E-$D, in which E denotes a projection operator defined 



where S{K) is a positive regularization parameter (not a (5-function!) and we 
have assumed that is self adjoint. As a final step, the parameter 

5{K) is reduced as much as required, and, in particular, when some second- 
class constraints are involved, 5{fi) ultimately remains strictly positive. This 
general procedure treats all constraints simultaneously and treats them all 
on an equal basis. 

C. Why ultralocal fields are relevant 

In quantizing any theory with constraints, including reparametrization in- 
variant theories, we invariably follow the rule: Quantize first, reduce second, 
also used by Dirac [^. [Remark: Some other quantization procedures re- 
duce (i.e., impose constraints) first and quantize second, a scheme that may 
lead to different and generally incorrect results, especially when the physical 
phase space (the quotient of the constraint surface by the gauge transforma- 
tions) is non- Euclidean 0. We do not comment further on these alternative 
procedures.] It must be strongly emphasized that in the initial quantization 
phase of this pair of operations, one must remain neutral toward, or even 
better, bhnd to any specifics of the constraints to be imposed. In particu- 
lar, in the initial, quantization phase, reparametrization invariance requires 
that basic fields are statistically independent at any spatial separation; cor- 
relations between fields that are spatially separated, and which originate 
from the particular physics of the theory under discussion, will arise only 
after the relevant constraints are imposed. Before the constraints are ever 
discussed, the nature of the field operators is ultralocal, the name given to 
field-operator representations which are, in fact, statistically independent for 
all distinct spatial points. The primary representation already presented in 
P-I for the basic local affine quantum field operators, and before constraints 
have been introduced, is for this very reason ultralocal in character. In the 
present paper, we shall discuss the relevant constraints for gravity and argue, 
to the extent possible at present, that the imposition of the constraints leads 



by 
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to field operator representations that are no longer ultralocal in character. 
Indeed, we shall argue that all traces of an initial ultralocal representation 
disappear when the constraints are fully enforced. [Remark: As an illus- 
tration of these general procedures, it is important at this point to observe 
that a relativistic free field of mass m (and, in addition, well-defined inter- 
acting fields) can be quantized starting from a reparametrization invariant 
formulation and with an initial field operator representation that is ultralocal 
in nature After imposing the appropriate constraint, it may be shown, 
by a suitable modification of the reproducing kernel — see Sec. IV for a brief 
discussion — how the conventional and nonultralocal field operator represen- 
tation for the relativistic free field of an arbitrary mass m emerges, as well 
as how the conventional propagator for such a system also arises.] 

D. Outline of remaining sections 

In Sec. II, we discuss in some detail the formulation of the initial stage of 
quantization in the absence of the gravitational constraints. In Sec. Ill, the 
gravitational constraints are introduced, and a novel, but generally familiar, 
functional integral representation of the desired expressions is developed. 
Finally, in Sec. IV, and aided by the use of several analogies, we present a 
lengthy discussion of the virtues, as we see them, of the specific functional 
integral representation developed in this article. 

II. REPRODUCING KERNEL FOR THE 
ORIGINAL HILBERT SPACE, AND ITS 
FUNCTIONAL INTEGRAL 
REPRESENTATION 

A. Basic field operator representation and the 
original Hilbert space 

For reasons briefly sketched above, we adopt as basic local field operators 
the positive-definite, local self-adjoint metric operators gab{x), and the lo- 
cal self-adjoint "scale" operators vr^(a;), which satisfy the affine commutation 
relations ([1|) in the "original" Hilbert space S). Here x = {x"}^^]^ denotes 
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the coordinates (each with the dimension of length) of a point in a classi- 
cal topological space S, which, for the sake of discussion, we may assume 
is topologically equivalent to M^. The term S{x,y) which enters the affine 
commutation relations is a Dirac 5-function "scalar density" with dimensions 
{length)~^ . Moreover, the affine commutation relations uniquely tell us that 
Tf^ix) transforms as a mixed- valence tensor density of weight one and has the 
engineering dimensions of an action (due to h) times {length)~^, i.e., M/LT 
in terms of mass (M), length (L), and time (T). We require that gab{x) 
transform as a covariant tensor of rank two, and that gab{x) is dimensionless. 

An auxiliary — but temporary — structure 

For purposes of the present section, we shall need to augment the classical 
topological space S with one additional structure, namely a volume form. 
Specifically, we adopt a real, positive scalar density of weight one, b{x), 
< b{x) < oo, X G iS, with which we may define, at each point, an invariant 
volume form dV = b{x) d^x. Additionally, in keeping with its transformation 
properties, we further insist that h{x) has the dimensions of L~^. 

The introduction of the auxiliary structure represented by b{x) is re- 
quired before the constraints are introduced, but whatever choice is made, it 
will disappear completely after the constraints are fully enforced. 

More explicitly, several arguments are offered below to justify the intro- 
duction of the auxiliary structure represented by the function b{x), x E S, 
prior to the introduction of the constraints. In Sees. Ill and IV, we will 
observe that we expect all the irrelevant freedom present in b{x) to be 
"squeezed out" by the constraints and replaced with another, nonarbitrary 
structure with the same dimensions. This is exactly the process observed 
in other cases that have been explicitly dealt with, and the particular case 
of a reparametrization invariant relativistic free field will be summarized in 
Sec. IV i. 

Original reproducing kernel 

Dual to the local operators gab{x) and ^^(x), we next introduce two real, 
c-number functions n°-^{x) [= n^°'{x)] and j^ix), which, initially, may be 
taken as smooth functions of compact support. Here tt""^ transforms as a 
contravariant tensor density of rank two and has dimensions M/LT, while 
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7c transforms as a mixed tensor and is dimensionless. With \ri) — called the 
fiducial vector — an as yet unspecified unit vector in ^, we consider the set of 
unit vectors (in units where ^ = 1) each of which is given by 

It,-, 7) = e*^ Tr''''ix)gab{x)(fx ^-i J -ti{x)nl(x)(f'x ^ 

As 7r"^ and 7^ range over the space of smooth functions of compact support, 
such vectors form a set of coherent states. 

The complex functional (vr", 7"|'7r', 7') formed by the inner product of 
two such coherent states will be a functional of fundamental importance 
in the present study of gravity. In particular, the functional (vr", 7"|7r', 7'), 
whatever form it takes, is manifestly a positive-definite functional that fulfills 
the defining condition that 



J 

E 

j,k=i 



a*ak (7rj,7j|7rfc,7fc) > (5) 



for general sets {aj} and {vTj, 7j} for any J < 00. Furthermore, (vr", 7"|7r', 7') 
is always a continuous functional in some natural functional topology, e.g., 
a topology defined by the particular expression itself 0. As a continu- 
ous, positive-definite functional, it follows that we may adopt the expres- 
sion (vr", 7"|7r', 7') as a reproducing kernel, and use it to define an associated 



reproducing kernel Hilbert space C |T^. Let two elements of a dense set of 
elements in C be given by 



K 

<P{t^,7) = ^(^k{Tr,l\7rk,7k) , K<oo, (6) 



k=l 



where {Tf^k,'lk}k=i denotes another independent set of (real) fields. These are 
continuous functionals of the fields vr and 7. As the inner product of these 
two elements we adopt 



J K 



(V^, = 2^ 2^ «j /9fc(7rj, 7jkfe, 7fc) • (7) 
j=i k=i 
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We complete the space of functions by including the limit point of all Cauchy 
sequences in the norm \\ip\\ = [ip^ipY^'^ . 

The result of the above construction is the (separable) reproducing kernel 
Hilbert space C composed of bounded, continuous functionals. Moreover, the 
Hilbert space C provides an especially useful functional representation of our 
original Hilbert space, which was referred to as the "primary container" in 



B. Choice of the fiducial vector and exphcit form of 
the reproducing kernel 

As argued in Sec. I, the representation of the basic field operators Qah and tt^ 
must be ultralocal prior to the introduction of any constraints. To fulfill this 
requirement, it is necessary that 



(7r",7'V',7') = exp{-/6(a;) AL[7r"(x),7"(a;);7r'(a;),7'(x)]} (8) 



for some dimensionless scalar function L. This function is determined by the 
representation of the affine field operators and the fiducial vector \ri), and as 
minimum conditions we require that 



Here, g{x) = {gab{x)} is a fixed, smooth, positive-definite metric function de- 
termined by the choice of \ri) (see Jl]]). Whether S is compact or noncompact, 
the choice of g{x) will determine the topology of the space-like surfaces under 
consideration; if S is noncompact, then g{x) also determines the asymptotic 
form of the space-like surfaces under consideration. 

For reasons to be offered below, we choose I?;) so that the overlap function 
of two coherent states is given (when ^ = 1) by 



P-I. 



{v\9ab{x)\r]) = (jabix) , 



(9) 
(10) 



(vr",7'V',7') 




X In 



det{\[g"^\x) + g'^'^x)] + \ib{x)-^[n'"'\x) - n"'\x)]} 
{det[^'"*^(x)] det[^'"^(x)]}V2 



) 



.(11) 
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Several comments about this basic expression are in order. 

Initially, regarding (|ll|), we observe that 7" and 7' do not appear in the 
explicit functional form given. In particular, the smooth matrix 7 has been 
replaced by the smooth matrix g which is defined at every point by 

g{x) ^ e^(-)/2 ~g{x) e^^-)"/^ ^ ^ (^g) 

where 7(2;)^ denotes the transpose of the matrix 7(x). Observe that the 
so-defined matrix {gab{x)} is manifestly positive definite for all x. The map 
7 — *■ 5^ is clearly many-to-one since 7 has nine independent variables at each 
point while g, which is symmetric, has only six. In view of this functional 
dependence we may denote the given functional in by {tt" , g" \7c' , g') , and 
henceforth we shall adopt this notation exclusively. 

Single afRne matrix degree of freedom 

An elementary example of the notational change from 7 to (7 may be seen 
quite directly in what occurs for a single affine matrix degree of freedom (in 
contrast to a field of such degrees of freedom). To that end, and following 
Ref. 11 closely, we introduce the Lie algebra for affine matrix self-adjoint 
operator degrees of freedom composed of the symmetric 3x3 matrix {a at} 
and the 3x3 matrix {/t^}, which together obey the affine commutation 
relations [cf., (|l])] 

[o-ab, f^d] = ^li^CTdb + S^CTad) , (13) 
[(7ab, (^cd] = . 

We choose the faithful, irreducible representation for which the operator 
matrix {aab} is symmetric and positive definite, and which is unique up 
to unitary equivalence. Furthermore, we choose a representation which di- 
agonahzes {aab} as k = {kab}, which we refer to as the ^-representation. 
In the associated representation space, and for arbitrary real matrices 
F = {F""^}, F'"^ = and B = {B^^}, it follows that 

U[F,B]ij{k) = e^^"'""^ e-'^>^{k) 

= {dei[S]fe'^'''^-'ilj{SkS^), (14) 
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where S = e ^/"^ = {S^} and (SkS'^)ab = S^kcdS^. The given transformation 
is unitary within the inner product defined by 

J ip{ky'4;{k)dk , (15) 

where dk = Ila<bdkab, and the "+" sign denotes an integration over only that 
part of the six-dimensional A;-space where the elements form a symmetric, 
positive-definite matrix, {kab} > 0. To define coherent states we choose an 
extremal weight vector, 

r]{k) = C(det[A;])^-^ ^-^^^[0-^1^] ^ ^^g^, 

where (3 > 0, G — {Gab} is a fixed positive-definite matrix, C is determined 
by normalization, and tr denotes the trace. This choice leads to the expec- 
tation values 

{'nWab\r]) ^ j r]{k)*kab'n{k)dk^ Gab , (17) 

Hk^M- j^r){kr>^ar){k)dk^Q. (18) 

In the /c-representation, it follows that the affine matrix coherent states are 
given by 

{k\F,B) = C(det[^])2(det[5A;5^])^-ie^*--[^'=le-''*--[^"^^^"l . (19) 

Observe that what really enters the functional argument is the positive- 
definite matrix G~^ = S'^G~^S where we set G = {Gab}- Thus without 
loss of generality wc can drop the label B (or equivalently 5") and replace it 
with G. Hence the affine matrix coherent states become 

{k\F, G) = C"(det[G'-^])''(det[A;])''-^ e'''^^''^ e"^*"^'^"'*^! , (20) 

where C — C(det[G])^ is a new normahzation constant. It is now straight- 
forward to determine that 

{F", G"\F\ G') = j (F", G"\k) {k\F\ G') dk 

{dei[G"-^]dei[G'-^]Yl'' 1^/3 
Aei{\[{G"-^ + G'-^) + i(5-\F" - F')]}\ ' ^ ' 

In arriving at this result, we have used normalization of the coherent states 
to eliminate the constant C. 
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Lattice construction 



Suppose now that we consider an independent lattice of such matrix degrees 
of freedom and build the corresponding coherent state overlap as the product 
of ones just like (0). Let n label a lattice site and let nGN, which in turn 
is a finite subset of 1? . In that case the coherent state overlap is given by 



(F",G"|F',G')n 



{det[Gf,^]det[G;„/]} 



nGN 



(22) 



As our next step we wish to take a limit in which the number of independent 
matrix degrees of freedom tends to infinity in such a way that not only does 
the lattice size diverge but also the lattice spacing tends to zero so that, 
loosely speaking, the lattice points approach the points of the space S. In 
order for the limit to be nonzero, it is necessary that the exponent /?[„] — >■ 
in a suitable way. To that end we set 

An]=^[n]A, (23) 

where A has the dimensions L^, and thus 6[n] has the dimensions L~^. In 
addition we need to let = vrj!^j A, and we rename Gab[n] as gabin] and call 
the matrix elements of Gj^J by g^^y With these changes (|2^) becomes 



{7r",g"\7T',g% 



{det[^?f;f]det[^?-^]}V2 



neN 



det{i[(^7ff + ^7S)+.6f4(7r-''-7r-f)]} 



2br„i A 



(24) 



Finally, we take the hmit as described above and the result is given by (pH]), 
namely. 



(tt ,g ,g} 



exp 



X In 



b{x) d?x 



/det{|[^"'^''(a;) + g"'\x)] + \ih{x)-^[T:"''\x) - n'^^x)]}^ 
V {det[^'"^^(x)] det[^'"''(x)]}i/2 , 



.(25) 
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In this way we see how the continuum resuh may be obtained as a hmit 
starting from a collection of independent affine matrix degrees of freedom. 
The necessity of ending with an integral over the space S has directly led to 
the requirement that we introduce the scalar density b{x). 

C. Additional arguments favoring b{x) 

As a further general comment about ( pT]) or (|25|) we observe that (vr", 5'"|vr', g') 
is invariant under general (smooth, invertible) coordinate transformations 
), and we say that the given expression characterizes a dijfeo- 
morphism invariant realization of the affine field operators. This property 
holds, in part, because b{x) transforms as a scalar density in both places that 
it appears. Thus b{x), which has the dimensions of L^^, plays an essential 
dimensional and transformational role in each place that it appears. Note 
that if h is explicitly introduced into ( pS]) or (0), it appears only in the 
change [k""'''{x) — 7r'"*(x)] [k""'''{x) — n'^'' (x)] / h. If one insisted on build- 
ing an acceptable ultralocal positive-definite functional using only 7r"°''(a;), 
g""'^{x)^ 7t'"-^{x), g'"'^{x), and h, then that construction would not appear to 
be possible simply on dimensional grounds. 

As another argument for the appearance of b{x), we observe that (^) 
involves not only the c-number fields vr and g (or 7), but the particular 
representation of the local operators gab and vr^ as well as the choice of the 
fiducial vector \ri). It is entirely natural that the function b{x) may emerge as 
a needed functional parameter in defining the operator representation and/or 
the vector \r]), and this property is explicitly illustrated in P-I. 

As a final argument for the appearance of the scalar density b{x), we 
briefiy recall properties of the local operator product for ultralocal affine 
field operators [|1|. In particular, the formal local product reads 

gab{x)gcd{x) = 6{x,x)Eabcd{x) + l.s.t. . (26) 

Here Eabcd{x) is a local fourth-order covariant tensor density operator of 
weight —1, 6{x,x) is a divergent multiplier with dimensions L~^, which arises 
when the "scalar density" delta function 6{x, y) is evaluated at coincident 
points, and l.s.t. denotes "less singular terms". Before adopting the proper 
local operator product, we introduce a scalar density b{x), < b{x) < 00, 
with dimensions L~^, and consider 

gab{x)gcd{x) = b{x)[b{xy^ 6{x,x)]Eabcd{x) + l.s.t. . (27) 
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Finally, we choose 



[gab{x) gcd{x)]R = b{x) Eabcd{x) {= gab{x)gcd{x)/[h{x) ^6{x,x)]) (28) 

as the proper renormalized (subscript R) local operator product. (Limits in- 
volving test function sequences offer a mathematically precise construction.) 
The given choice leads to a local operator that transforms as a tensor in the 
natural fashion and, moreover, carries the natural engineering dimensions. 
To achieve this desirable property in a local product requires the introduc- 
tion of an auxiliary scalar density h{x). 

For additional properties regarding local products of the relevant affine 
field operators, see P-I. Essentially, all these properties are almost entirely 
based on the analysis of local operator products in scalar ultralocal field 



theories, an analysis which is described in detail, for example, in [|12 |. 

We have offered several reasons for the appearance of the scalar density 
h{x) at the present stage of the analysis. However, we emphasize once again 
that h{x) will disappear when the constraints are fully enforced, whatever 
choice was originally made. An example of the process by which this funda- 
mental transformation takes place is presented in Sec. IV. 

D. Functional integral representation for the 
coherent state overlap functional 

For further analysis, especially when we take up the issue of introducing the 
constraints in Sec. Ill, it is useful to introduce an alternative representation 



of the functional expression (^51) . The alternative representation we have in 
mind is that of a specific functional integral, which, indeed, has already been 
introduced in P-I. That such a representation should exist is an immediate 
consequence of the fact that ( pSj ) fulfills a complex polarization condition, 
which then leads to (p5D being annihilated by a negative, second-order func- 
tional derivative operator. Exponentiating this operator, times a parameter 
z/ > 0, letting the resultant operator act on general functionals of vr" and g" , 
and including any necessary i/-dependent prefactor, will, in the limit v oo, 
lead to a dense set of functions in the reproducing kernel Hilbert space C. 
Alternatively, letting the same operator act on a suitable (5-functional will 
lead to the expression (pSf). The functional integral of interest arises in this 



last expression by introducing the analog of a Feynman-Kac-Stratonovich 
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representation. The mathematics behind these foregoing several sentences 
is well illustrated in P-I for both a simple, single affine degree-of- freedom 
example as well as for the affine field theory. 

The result of the operations outlined above leads to a functional integral 
representation for {ir" , g"\'n'' , g') in (^51), which is given (for h = 1) hj 



I II II \ I i\ 
(tt ,g |vr 



exp 
X In 



h{x) (fx 



/det{i[^""^(a;) + g'^^x)] + ^ib{xy^[Tr"''''{x) - Tr"'\x)]}\ 
V {det[g"^^{x)] det[^"**(a;)]}V2 ) 

= lim 77 ,y / expl-iJgabTr"'' d^xdt] 

X exp{~{l/2u)J[b{xr'g,b9cdTr'''7T''^ + h{x)g''' g'^g^cgda] d'x dt} 
X n,,j n,<f, dTT''\x, t) dgabix, t) . (29) 

Here, because of the way the new independent variable t appears on the right- 
hand side of this expression, it is natural to interpret t, 0<t<T,T>0 
as coordinate "time" . The fields on the right-hand side all depend on space 
and time, i.e., gab = gab{x,t), gab = dgab{x,t)/dt, etc., and, importantly, the 
integration domain of the formal measure is strictly limited to the domain 
where {gab{x, t)} is a positive-definite matrix for all x and t. For the boundary 
conditions, we have 7r'"^(x) = 7r'**(x, 0), g'ab{x) = gab{x, 0), as well as 7r'"*^(x) = 
7r"*(x, T), g'abix) = gab{x, T) for all x. Observe that the right-hand side holds 
for any T, < T < oo, while the middle term is independent ofT altogether. 

Although the functional integral on the right-hand side is formal it nev- 
ertheless conveys a great deal of information. Let us first examine it from 
a dimensional and transformational standpoint. As presented, h = 1; to see 
where h would appear we may simply replace each tt"^ by n"'^/ h. With u hav- 
ing the dimensions of and AT,^ absorbing any remaining dimensions from 
the formal measure, then the right-hand side of (^91) is dimensionally satis- 
factory. From the point of view of (formal) transformations under coordinate 
changes, it is clear, with AT,^ transforming appropriately, that the right-hand 
side is formally invariant under coordinate transformations involving the spa- 
tial coordinates alone. (Remark: A discussion about transformations of the 
right-hand side under spatially dependent transformations of the time coor- 
dinate has been given in P-I and is not repeated here. It is clear that the 
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result of the limit z/ — > oo on the right-hand side must be invariant under all 
such transformations simply because the middle term is independent of the 
time variable altogether!) 

As presented — and indeed as originally derived — the result of the func- 
tional integral (the middle term) came before the functional integral repre- 
sentation of that result (the right-hand side). However, we can also interpret 
( |29D in the opposite order, that is, to presume that the functional integral 
(right-hand side) is primary and that the answer (the middle term) is the 
result of evaluating the functional integral. This perspective encourages us 
to examine the expression in the integrand of the functional integral some- 
what more carefully from a traditional standpoint. We first observe that the 
formal, flat part of the measure has the expected appearance of the canoni- 
cal measure for a conventional, canonical functional integral quantization of 
gravity. The phase factor contains an acceptable classical symplectic poten- 
tial term in a formal functional integral for gravity in which the rest of the 
classical action — the terms involving the constraints — are absent. (Remark: 
This characterization is, of course, quite appropriate since we still are in the 
first phase of our dual approach: quantize first, reduce second. The patient 
reader will be rewarded with the addition of the expected constraints and 
Lagrange multiplier terms in the next section.) The second, zz-dependent fac- 
tor in the integrand serves as a regularizing factor for the functional integral. 
Formally, as — > oo, such a factor disappears from the integrand leaving the 
expected (pre-constraint) formal functional integral integrand, such as it is. 
However, the zz-dependent term plays a fundamentally important role within 
the integral itself since it literally serves to define the functional integral. 

It is important to make this last point quite clear, and for that purpose 
we temporarily discuss the formal expression (with h = 1 again) 



dfi''{n,g) 

= exp{-(l/2zz) /[6(a;)"'^7a6^7cd^'^7r'^» + h{x)g''''g'^''g,,ga^\ d\ dt} 



We assert that for fixed b{x), < b{x) < oo and fixed u, < u < oo, this ex- 
pression characterizes a bona fide, countably- additive, positive measure, fi", on 
the space of generalized functions vr"* = 7r°''(a;,t) and gab = gab{x,t), where, 
for any nonvanishing and any nonvanishing, nonnegative test function 
f{x) > 0, the positive-definite matrix condition J f{x)u"' gab{x,t)u'' d^ > 




(30) 
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holds for (almost) all t, < t < T. The fields ir"'^ and gab satisfy the boundary 
conditions at t = and t = T given previously. The factor J^i, is adjusted 
so that the measures /i'^ form a semi-group with respect to combining time 
intervals, e.g., — T, T > 0, and then T T + T' , T' > 0, being equivalent 
to ^ T + T'. If {hp{x)}'^^ denotes an orthonormal set of test functions 
defined so that 

Jhp{x)hq{x)b{x) d\ = 5pg , (31) 
b{x)E^^,hp{x)hp{y)=6{x,y), (32) 

then we assert that there exist finite, nonzero constants Np for all u and all 
P e {1,2,3,...} such that 

N^p j eM-^K=i!gaKp)mt){t) dt] dfi''{7r,g) (33) 

is well defined. In this expression, 

gab{p){t) = Jhp{x)gabix,t)b{x) dh , (34) 
7r^p^)(t) = Jhp{x)Tr''\x, t) d^x . (35) 

Moreover, the set of constants {Np} may be chosen so that 

(tt ,g \7T ,g} 

= lim lim N^ I exp[-iE^^Jgab(p)it)n?\{t) dt]diJ,''{7r,g) . (36) 

P^oo i/—*oo J 

This is one of the ways that the formal functional integral (|29D can be given 
a rigorous meaning. 

There is another way to give rigorous meaning to (^) that we would 
also like to discuss. In this procedure we use a spatial lattice but keep the 
time variable t continuous. This regularization scheme takes us back to the 
lattice construction given earlier [cf., (p^j, except now we add a phase-space 
path integral representation as well. For present purposes we again introduce 
the symbol A, with dimensions L^, to denote a uniformly (coordinate) sized, 
small spatial cell. Then the lattice regularized path integral expression for 
{7i",g"\Tr',g')f^ in (E3) is given by 



X n„,t n,<, rfTrg (t) dgab[n] (t) , (37) 
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where n G N, which itself is a finite subset of 71?. Here ttj^j, gah\ri\^ and 
6[n] represent average field values in the cell n, and h^^^ = l/&[n]- Next, we 
observe that there is a countably-additive, pinned, Brownian-motion measure 
formally defined by 

= Af^ exp{-^i:nJ[bl^]9ab[n]9cd[n]Tr'(^]^l^ + &[n]^g^[n]^6c[n]^da[n]] A rft} 

X n„,t Ua<b rfvrg (t) dgab[n] (t) (38) 
SO that (0) becomes 

iV^ ! e^'''-^'^'^-^<'"'Uii-^{n,g) , (39) 



where {N^ } is a set of finite constants. Moreover, these constants may be 
chosen so that 

(7r",/|7r',(7')= lim \im f e-'""-^ '^'^-^^'^''U^UT^^g) . (40) 

N^oo v-*oo J 

Here, in the last step, the limit N — oo means that A ^ and N — 
in such a way that all points x G 5 are reached in a natural way. Observe 
that the presence of the continuous-time regularization factor in the formal 
functional integral for the entire space has controlled the spatial lattice reg- 
ularization in a clear and natural fashion; although possible to introduce, no 
temporal lattice has been required to obtain a well-defined expression. 

The present usage of well-defined, phase-space measures to define func- 
tionals integrals in the limit that the "diffusion constant" parameter (z/) di- 
verges is part of the general program of continuous-time regularization [|l^ . 
It is noteworthy that the use of a suitable phase-space metric to control the 
Brownian motion paths invariably leads to a coherent-state representation for 
the resultant quantum amplitude. These remarks conclude our brief excur- 
sion into a rigorous discussion of the functional integrals of present interest. 

Another prospective regularization 

Let us return to the formal functional integral ( ^9]) and examine that expres- 
sion with regard to the scalar density h{x). Superficially, in the limit v ^ oo 
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in which the regularizing term proportional to {l/2v) in the exponent of the 
integrand formally vanishes, one might naively expect that the result of the 
integral would be independent of the function h{x). But, no, that naive ex- 
pectation is false since the result of that integral and subsequent limit, i.e., 
the middle term of (p9[) , evidently depends importantly on h{x). At first 
glance, that dependence seems highly unnatural. It may be thought — as 
the author did for a number of years [0 — that an alternative regularization 
expression may be more "natural" and would therefore be preferable. In 
particular, it was thought that the expression 

exp{-(l/2z/)/[<7-i/2(7,,^?,,7r^%'^'^ + g^'^g'^'g^^g^cQda] d\ dt] , (41) 

where g = det[gab], and which involves a different phase-space metric, was 
"better". However, in light of the discussion in the present paper, it is now 
evident that this expression is not even dimensionally consistent! This defect 
could be rectified by the introduction of a positive constant, which we may 
call b, that stands in the place of the present b (next to g^^"^) and carries the 
dimensions L~^. Indeed, we could also introduce in place of 6 a positive scalar 
function b{x) with dimensions L~^. In any case, the failure of (^) purely on 
dimensional grounds, and the necessity thereby of introducing some sort of 



auxiliary dimensioned parameter (or function) surely renders (¥1) far less 
"natural" than had been previously assumed. 

As another possible argument against (PT|), we note that if that form of 
a proposed regularization was used in an expression like (^) in place of the 
present form of regularization, there is absolutely no guarantee that the re- 
sult will describe a reproducing kernel with other than a one- dimensional 
reproducing kernel Hilbert space, which is the general result for a "random" 
choice of phase-space metric. The fact that the present form of (|29| ) generates 
a suitable, infinite-dimensional reproducing kernel Hilbert space is a funda- 
mentally important feature, which, in our case, is a consequence of having 
started with an appropriate reproducing kernel to begin with. 

Although the author has not foreclosed any possible interest in a 6 mod- 
ified version of the regularization (|41|), all present indications point to the 
version (|30D that features the scalar density b{x). This shift of allegiance has 
also been bolstered by the realization that the role of b{x) is confined to the 
initial phase of quantization and that b{x) will disappear entirely after the 
constraints are fully enforced. On the strength of this argument, it is the 
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version based on ( |30D and not (^Tj) that is analyzed in the remainder of the 
present paper. 



III. INTRODUCTION OF THE 

GRAVITATIONAL CONSTRAINTS 

A. Key principles in heuristic form 

There are four gravitational constraint functions, the three diffeomorphism 
constraints 

Haix) = -2nl^,{x) , (42) 

where "|" denotes covariant differentiation with respect to the spatial metric 
gab, and the Hamiltonian constraint, which, in suitable units (i.e., c^/G = 
IGtt), reads 

H{x) = ^7(x)-i/2[<(x)7r^(x) - i<(a;)<(a;)] + ^^(x)V2(3)/2(^) , (43) 

where g{x) = det[gab{x)] and ^^^R{x) denotes the scalar curvature derived 
from the spatial metric [|T6|. Classically, these constraint functions vanish, 
and the region in phase space on which they vanish is called the constraint 
hypersurface. 

It is instructive to evaluate the classical Poisson brackets between the 
constraint fields. For this purpose, we enlist only the basic nonvanishing 
Poisson bracket given by 

{gabix), n^\y)} = '^{6^Jt + ^X) K^, v) ■ (44) 

It follows that 

{Ha{x),Hb{y)} = 5^a{x,y)Hb{x) - 5^b{x,y)Ha{x) , (45) 

{Ha{x),H{y)} = 5,a{x,y)H{x), (46) 

{H{x),H{y)} = 5^a{x,y)g^\x)Hb{x) . (47) 

In these expressions, 5^a{x-,y) = d6{x,y)/dx'^, which transforms as a "vec- 
tor density". It is clear that the Poisson brackets of the constraints vanish 
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on the constraint hypersurface because the right-hand sides of (^)-(^) all 
vanish there, i.e., when Ha{x) = = H{x) for all x G iS. This vanishing 
property of the Poisson brackets is characteristic of first-class constraints. 
Often, the Poisson bracket structure of constraints is that of a Lie bracket in 
which case such constraints are referred to as closed first-class constraints. 
In order for the Poisson bracket structure to be that of a Lie bracket, it is 
necessary that the coefficients of the constraints on the right-hand side in- 
volve no dynamical variables. Taken by themselves, we note from (^Sj) that 
the three diffeomorphism constraint functions form a set of closed first-class 
constraints. However, because of the last equation, p7|), it is clear that the 
complete set of four gravitational constraint functions do not have a Poisson 
structure which is that of a Lie algebra, and consequently the gravitational 
constraints are said to form an open first-class system of constraints. Such a 
situation does not automatically imply trouble in the corresponding quantum 
theory, but significant difficulties do arise in a number of cases. Quantum 
gravity is one of those cases. 

Let us proceed formally in order to see the essence of the problem. Sup- 
pose that Ha{x) and T-C{x) represent local self-adjoint constraint operators 
for the gravitational field. Standard calculations lead to the commutation 
relations (with ^ = 1) 

[na{x),nb{y)] = I [5 a(x, y) n^[x) - y) Uaix)] , (48) 
\na{x),n{y)\=%b^a{x,y)n{x) , (49) 
[7^(x),7^(y)] = zi5,(a;,y)[^'^^(a;)Hb(a;)+7^,(x)^"^(a;)] , (50) 

where to ensure the Hermitian character we have symmetrized the right-hand 
side of the last expression. In the usual Dirac approach to constraints alluded 
to in Sec. I, one asks that \'^)phys = for all constraints. If we assert that 
'Ha{x)\ip)phys = and 'H{x)\ip)phys = 0, then consistency holds for the first 
two sets of constraint commutators, but not for the third commutator in 
virtue of the fact that it is almost surely the case that g°''^{x) \'ip)phys ^ ^phys, 
even if it were smeared. The expected behavior is somewhat like that of 
the single degree-of-freedom example where Q\ip) = and P\il)) = imply 
for standard Heisenberg operators that [Q,P] li') = il'ip) = 0, i.e., there are 
no nonvanishing solutions. This behavior is characteristic of second-class 
constraints, and as a consequence of our discussion we are led to conclude, 
from a quantum mechanical standpoint, that part of the gravitational con- 
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straints are second-class constraints |T^. For the projection operator method 
of constrained system quantization, however, second-class constraints cause 
no special difficulty and, in particular, they are treated in just the same 
way as first-class constraints, as already noted in Sec. I. (Remark: Some 
researchers prefer to modify the theory so as to eliminate the second-class 
nature of the gravitational constraints. Instead, we accept the second-class 
constraints for what they are.) 

Assuming that the constraint fields Ha{x) and H{x) are local self-adjoint 
operators, we could — as one of several different alternatives — proceed as fol- 
lows. Initially, besides the real, orthonormal set of test functions {hp{x)} 
introduced in Sec. II, let us introduce an additional set of real test functions 
{/p^(a:)}, A G {1,2,3}, with the following properties: 

6(a:)E- = r\x)6{x,y) , (51) 

/ fpA i^) fqBa{x)b{x)d^X = 6pg 6 AS , (52) 

where fqBa{x) = gab{x) fgj^{x). Regarding the set of functions {/^^(a;)}, the 
index A is a dreibein index, while the index a is a three-space vector index. 
With the help of these sets of expansion functions, let us introduce 

Hip) A = JfpAix)na{x) d\ , (53) 

'H(p) = Jhp{x)n{x)d^x , (54) 

each for 1 < p < oo, and similarly for other vector and scalar functions. 
In this form as well, part of the constraint operators are second class. To 
accomodate all these constraints we introduce a set of projection operators 
defined for all P G {1, 2, 3, . . .} and given by 

Ep = E(X|, < 6{hf) , (55) 

x| = s;^, {n^p)Ar + nfp) ] . (56) 

As defined, the projection operators Ep are regularized and they serve to 
define regularized physical Hilbert spaces S^phys = Epi). These regularized 
physical Hilbert spaces may, in turn, be characterized by their own repro- 
ducing kernels 

(7r",/|Ep|vr',^'), (57) 
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and the regularized physical Hilbert spaces themselves may, therefore, be 
represented by the associated reproducing kernel Hilbert spaces. 

The final step in the present construction procedure would involve suitable 
limits to remove the regularizations. More familiar procedures to enforce the 
constraints are discussed in Sec. IV. 



B. Functional integral representation for the 
relevant projection operators 

In an earlier work , we have presented a very general procedure to construct 
the projection operator 1E(Sq,$^ < S{h)'^) by means of a universal functional 
integral procedure. In particular, it follows that 

]E(E„$2 < s{hf) = Jt e-^/^"W*" VR{\) , (58) 

where T denotes the time-ordering operator, {A°(t)}^^]^, < t < T, denotes 
a set of c-number "Lagrange multiplier" functions, and VR{X) denotes a 
formal measure on such functions. A suitable measure R may be determined 
as follows: First, introduce a Gaussian integral over the set {A"(t)} so that 

X-y y"Te~*i"^"(*^*"'^*e('/^^^/^"^" PA = e"'^^"*' . (59) 

Second, and last, integrate over 7 according to the rule 

r sinUSihf + C]} ^ , ^ . 

C^o+ L^oo 7r7 

The inclusion of the variable ( and the limit C — > 0+ ensures that we include 
the equality sign in the argument of E. Observe that this construction is 
entirely independent of the nature of the set of constraints 

Next, we continue to proceed formally in order to envisage how the pro- 
jection operator method may be used in the case of gravity. Adopting the 
foregoing analysis, we suggest that 

Ep = E(E^.i 2-' I EL, (W(rt a)' + H(p, I) 
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for an appropriately defined formal measure R. Observe that the integral 
over the measure R may well integrate over degrees of freedom that are not 
present in the time-ordered product (such as for p > P); however, there is 
no harm in doing so since R is already defined so that J VR = 1. We may 
also go one step further and assert that as P — oo we obtain the formal 
expression 

E = y Te-'/[^'"«»+^«]^'"''*I5i?(iV^iV) , (62) 

which, heuristically at least, realizes the projection operator that projects 
the original Hilbert space onto a correspondingly regularized physical Hilbert 
space. 



C. Functional integral representation of the 

reproducing kernel for the physical Hilbert space 

In Sec. II, we presented in ( P^D a continuous-time regularized functional in- 
tegral representation of the reproducing kernel (ir" , g"\7r' , g') for the original 
Hilbert space. The reproducing kernel for the (regularized) physical Hilbert 
space is given, in turn, by the expression {7[",g"\'E\n',g'). In order to give 
this latter expression a functional integral representation we first regard 

J[N''na + Nn]dh (63) 

as a time-dependent "Hamiltonian" for some fictitious theory, in which case 

= lim AT^ / exp{-i JlgabT^'"^ + N''Ha + N H] d\ dt} 

X exp{-{l/2u) J[b{x)-'gabgcdri'''7f'''' + h{x)g^''g''^g,,gda] d\ dt} 
X H,,i Ua<b dn''\x, t) dgabix, t) . (64) 

In this expression, there appear symbols Ha{7T,g) and H{7T,g) correspond- 
ing to the quantum operators Ha and H. Superficially, these symbols may 
(formally) be identified with the classical diffeomorphism and Hamiltonian 
constraint functions, in which case the expression (|6^) contains in its phase, 
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and up to a surface term, the full Einstein action. Thus (Q) comes ever 



closer to looking like a more traditional functional integral for gravity. 

However, before integrating over the functions A^"^ and and completing 
the story, we need to caution the reader that h is not zero (but rather one) and 
therefore the symbols Ha and H may not coincide with their usual classical 
expressions. All we can say at present is that Ha is a symbol for the operator 
Ha, a = 1,2, 3, and that if is a symbol for the operator TC. The connection 
between symbol and operator is implicitly contained in (|64D, and since, for 
the moment, the functions Na and A^ are general functions within our control, 
we may use that fact to assert that 

{n",g"\J[M%y)na{y) + M{y)H{y)]d'y\n',g') 

= limAT, / e-'I^'^"''"'''^''' f[M-{y)Ha{y,s) + M{y)H{y,s)]d]j 

X exp{-(l/2z/) j[h{x)-^gab9cd'f^'"'h^'' + Kx)g'^' g^'^g.^g^a] dhdt} 
xIl^^t'n^a<bd7T''\x,t)dgab{x,t) (65) 

for any smooth (test) functions M*^ and M and for any time s, < s < T. 
Equation (|65|) gives an implicit connection between symbol and operator 
for the present theory. We observe that the more traditional connection 
between symbol and operator that normally holds for Wiener-regularized 



coherent-state path integrals is unavailable in the present case since we 
are dealing with so-called weak coherent states for which no resolution of 
unity exists; see [|l], 0]. 

In addition, thanks to analyticity in the present case, the diagonal matrix 
elements of an operator uniquely determine the operator, and so we can 
also assert the connection between symbol and operator (given by setting 
n",g" = n',g'), 

J[M%y){7,\g'\naiy)\n\g') + Miy){n\g'\niy)\n\g')]d'y 

= limAT, / e'^f'^'^"'''"''''^''' f[M-{y)Ha{y,s) + M{y)Hiy,s)]d]j 

X e^v{-{l/2v)j[h{x)-^gabgcdT^''7^'"' + h{x)g'''g'''g,,gda\ d\dt} 

d'K''\x,t)dgab{x,t) , (66) 

again for smooth functions and M and any s, < s < T. We note that 
{t:' , g'\l-La{y)\T^' , g') and {-k' , g'\l-L{ii)\Ti' , g') denote still other symbols that are 
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often associated with the local operators Ti.a{x) and 7^(x), respectively. In 
(|66D , the notation f Qabi^"''^ dt means that only closed paths in phase space 
enter, i.e., just those paths for which the functions 

7r"*(x, 0) = 7r"^(a;, T) = 7r'"^(a;) , (67) 
gab{x, 0) = gab{x, T) = g'^^{x) (68) 

for all X E S. Note that a closed line integral in phase space involves just 
the symplectic form, and the result of the integral § Qab'^"''' d^x dt is invariant 
under any (smooth) change of canonical coordinates. 



Reproducing kernel for the physical Hilbert space 

We now complete the story by interpreting the otherwise arbitrary c-number 
functions N"- and N as Lagrange multiplier functions and integrating them 
out of (0). Since A^" and are not dynamical variables that must enter 
the formal phase-space functional integral measure in a prescribed way (i.e., 
as '^dpdq"), we are free to integrate them as we choose — and we choose to 
integrate them in such a way as to enforce the quantum constraints, at least 
in a regulated fashion. As explained above, one natural way to achieve our 
goal involves the formal integration measure 'DR{N°' ,N). 

Combining several steps previously described, we now assert that the 
reproducing kernel for the regularized physical Hilbert space has the phase- 
space functional integral representation given by 

{n\g"\JLW,g') 

lim AT,, I (,-ij[9ahi^''''+N''Ha+NH]cfxdt 
v^oo " J 

X exp{-(l/2z/) j[h{x)-^gabgcdi^''Tr'"' + h{x)g''^g''^g,,g,a] d'xdt} 
X n,,i Ua<b d7i''\x, t) dgabix, t)] VRiN'', N) . (69) 

In this final expression we have reached our primary goal, at least from a 
formal perspective. Despite the general appearance of (|69D, we emphasize 
once again that this representation has been based on the affine commuta- 
tion relations and not on any canonical commutation relations! Later, we 
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shall discuss a more careful definition of this formal expression along lines 
introduced in Sec. II, but for now let us examine (^) for its own sake. 

We first comment on the range of integration for the "lapse" variable 



N{x,t), a subject of recurrent interest |T^. Our view is that the range 
— oo < N{x,t) < oo is the proper range when quantizing the theory. After 
all, in the Hamiltonian viewpoint, space-time is a derived structure of the 
classical theory. In principle, the issue here is no more complicated than 
for the reparametrized one- dimensional free particle. For this example, the 
original classical action is / = J[pq — ^p'^]dt, where q = dq/dt, and has 
solutions p{t) = Po and q{t) = Pot + qo- The reparametrized version is 
given by /' = j[pq* + st* - \{s + \p'^)]dT, where q* = dq/dr, t* = dt/dr, 
and has solutions p(r) = po, s(r) = So = —^pI, tij) = \{a)da, and 
g(r) = t(r)po + go. The function A(r) is essentially arbitrary. If, for example, 
A(r) = 3r^ — 1, i.e., t(r) = r(r^ — 1), the solution seems "to go backward 
in time", but that interpretation gives to the variable r an unwarranted 
physical significance. The given solution is not wrong, it just repeats itself 
for a while. We can avoid a repeating behavior, e.g., by simply dropping 
the interval —1 < r < 1. No such issues occur if we require that A(r) > 
for all r, in which case r does indeed merit the name of "reparametrized 
time" . By analogy, the function N{x, t) — which we have loosely called the 
lapse function — only deserves that name when, in the classical solution space, 
we insist that N{x,t) > 0; otherwise it is just another Lagrange multiplier 
function, no more and no less. 

When one starts from a classical perspective, with its focus on physically 
relevant functions N{x,t) which are strictly positive, it is a conceptual leap 
to change to functions N{x,t) that can take on both signs |jl8|. However, 
when one starts from the quantum theory, as we have done, there is no such 
leap to make. 

As a second topic regarding (|69D we focus on its general structure. With 
Ha and H formally equal to the constraint functions of gravity (possibly up 
to terms in h), the action appearing in the phase factor is indeed appro- 



priate to gravity Moreover, the domain of integration is restricted to 
positive-definite metrics {gab{x,t)} > 0. Indeed, the particular i/-dependent 
regularizing phase-space metric in (|69D prevents the metric variable {gab} 



from escaping the positive-definite domain. 
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As a useful analogy, we note that the two-dimensional phase space metric 

(3-\^dp^ + (3q-^dq^ (70) 

is geodesically complete in the half-space {p, g) G M x ]R+, and when it is part 
of a Brownian motion measure, as in the formal expression [cf., (|38D] 

it automatically restricts the Brownian motion trajectories to the half-space 
R X M+. 

Lastly, we comment on the formal functional integration measure in (|69|), 
specifically for the Lagrange multiplier functions A^" and A^. The formal 
measure VR has been defined earlier and is unlike conventional measures 
chosen for such variables. As emphasized here, and elsewhere [Q, the mea- 
sure VR is designed to implement the quantum constraints — as befits a quan- 
tum theory — and it has not been selected to enforce the classical constraints. 
Observe well that we have not "blindly postulated" the functional integral 
(|69D, but instead it has been derived as a specific functional integral repre- 
sentation of the well-chosen quantum matrix elements on the left-hand side. 
Unfamiliar as the measure VR may appear, we maintain that VR — or some 
other measure equivalent to it — is the proper measure to choose to achieve 
our goal. Whether different treatments of the Lagrange multiplier functions 
that have been adopted by other workers are indeed equivalent or not to the 
use of VR is an interesting question, but it is not one we pursue here. 



IV. DISCUSSION 

In the preceding analysis, we have been strongly guided by the operator 
structure of the assumed theory of affine quantum gravity, and this discussion 
has led us to the formal functional integral representation ( ]69| ) for the desired 
matrix elements. As was previously the case [cf., (p9D], we now wish to turn 
( pUD around and adopt the formal functional integral as our starting point 
and, in effect, use that expression to define the reproducing kernel for the 
regularized physical Hilbert space. Specifically, for that purpose, let us focus 
on the formal expression 

hm 77^ [ e'^f [■^""'aab-N'^Ha-NH] d^xdt 

1/— »oo / 
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X exp{-(l/2i/) J[b{x)-'gabgcd7r''7r'''' + h{x)g''^g''^g,,gda] d'xdt} 
xVTrVgVR{N'',N) . (72) 

Note the change of the kinematic term, which simply amounts to a phase 
factor in the definition of the coherent states. Also we have introduced the 
common shorthand VnVg for the bracketed term in (|5UD. For the sake of 
discussion, we shall refer to (|72D as the "nonstandard expression" . Our goal 
in this section is to discuss the nonstandard expression and see what steps 
are necessary to give it a proper meaning. We shall do so in a three step 
procedure: First, we compare the "standard" (see below) and "nonstandard" 
expressions. Second, we discuss a regularization and its removal that ulti- 
mately involves the elimination of the scalar density b{x). Third, we examine 
the aspect of the problem that normally accounts for the perturbative non- 
renormalizability of the gravitational field. 



A. First look at the Nonstandard Expression 



It is interesting to compare (|72]) with what we refer to as the standard ex- 
pression for a phase-space functional integral for gravity. By the "standard 
expression" we mean the formal functional integral 

M y"e*/[-"'^"''-^"^"-^^l'^'"'^*p7rP(7PiV, (73) 

where VN = Ilx tdN{x,t) IladN'^{x,t). In several important ways, the stan- 
dard expression is very different than the nonstandard expression. Let us 
first comment on some of those differences. Much as it would be nice to 
think otherwise, one must recognize that the standard expression ( [73D is to- 
tally undefined as it stands; it is little more than a fancy way of writing 
X (oo). It begs for a definition as the limit of meaningful expressions [much 
as X (oo) may, for example, be defined as lim^-^o ^ x i'^/^) = but what 
set of meaningful expressions should be chosen in the gravitational case is 
far from clear. A lattice limit? But then, what form should the regularized 
lattice expressions take? Symmetry and covariance offer only limited guid- 
ance. In point of fact, this same question faces any standard phase-space 
path integral, even that for a single degree of freedom in which a conven- 



tional lattice definition — as originally envisaged by Tobocman |]19| — makes 



certain assumptions about the nature of the phase-space coordinates which 
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may or may not be true. [Remark: The skeptical reader is urged to pro- 
pose a lattice prescription to quantize the nonrelativistic free particle of unit 
mass by a phase-space path integral whose Hamiltonian H is expressed in 
canonical coordinates p and q such that H = \{p^ + q^)\ 

A further complication of the standard functional integral expression (ff^) 
arises from the unbounded nature of the formal integral / ■ ■ ■ T^N over the 
Lagrange multiplier variables. This choice of integration measure, which is 
designed to enforce the classical constraints and thereby reduce the classi- 
cal phase space before quantization, necessitates gauge fixing (to eliminate 
concomitant divergences) which reduces the classical phase space further to 
the physical phase space (at least locally) where each point labels a phys- 
ically distinct state. Quantization on the reduced phase space is formally 
aided by the introduction of some additional factor (e.g., a Faddeev- Popov 
determinant, or its analogue), which may well lead to significant (Gribov) 
ambiguities which require a substantial modification of the functional integral 
[pO| , 1^ , and give rise to serious problems (such as unitary violation) within a 



BRST formulation 121 



Nonstandard expression 



Let us raise similar issues regarding the nonstandard expression (72). Al- 
though ( |72D is formal as it stands, it can, to a considerable extent, already 
be regarded as "nearly" well defined. As noted in Sec. II, we can combine sev- 
eral factors together to make, for all finite a positive, countably-additive, 
pinned measure on generalized functions. What makes (^) not well de- 
fined is the fact that the formal integrand does not constitute an integrable 
function with respect to that measure. We have already encountered that 
problem in Sec. II before any constraints were introduced, and we found that 
we could overcome that problem by regularizing the integrand and removing 
that regularization as a final step. Superficially, the same property holds 
when the constraint functions are present (say for fixed Lagrange multiplier 
values), save for one very important distinction (involving the field opera- 
tor representation) which we shall address in the second point of discussion 
below. 

Regarding the integration over the Lagrange multiplier variables, we em- 
phasize the vast difference afforded by the projection operator method. First 
and foremost is the fact that no gauge fixing is introduced, no ghosts are 
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used, no Faddeev-Popov determinant (or its analogue) arises, and conse- 
quently, no Gribov ambiguities can exist. These properties arise, largely, 
because j'DR{N°',N) = 1, while jT>N = oo. The difference here could 
not be greater, and it arises because in the former case one quantizes first 
and reduces second, while in the latter case one reduces first and quantizes 
second. Except in basically trivial cases, the second option is fraught with 
substantial obstacles 0. 

One of the most significant differences between the standard and the non- 
standard expressions refers to the representation of the quantum mechanical 
amplitudes that is involved. For the standard expression, it is usually as- 
sumed that (^) leads to a representation in which the metric field operator 
gab{x) is diagonalized and thus is sharply represented. In combination with 
any needed auxiliary factor in the functional integral, diffeomorphism in- 
variance suggests that ([7^) depends only on the "geometry" of the initial 
and final three-surfaces, and not on the details of any specific metric ex- 
pressions jT^ , |T8[] . An analogous view is also prominent in the associated 
"loop quantum gravity" in which, e.g., bras and kets depend only on knot 
invariants as labels of "physically" distinguishable states p2[. In contrast. 



the representation afforded by the nonstandard expression (^) is that of a 
coherent- state representation, which depends on smooth metric gab and mo- 
mentum 71°'^ fields, that represent not sharp operator (eigen)values but mean 
operator values, e.g., 

{n,g\n^\x)\n,g) = n''\x), (74) 
{T^,g\gabix)\n,g) = gabix) , (75) 
{7r,g\7rl{x)\7r,g)=gUx)7i'\x) . (76) 

Note well, that besides the local self-adjoint metric gab{x) and scale field 
^^(x), we have used the local symmetric (but not self-adjoint) momentum 
operator 7r"^(x) in the first of these expressions. These expectations are not 
gauge invariant, nor should they be, since they are taken in the "original" 
Hilbert space where the constraints are not fulfilled. The gauge invariant 
part of the metric field, for example, and in so far as the regularized physical 
Hilbert space is concerned, is determined by the matrix elements 

{Ti'\g"\I^gab{x)'EW,g') , (77) 

which is an expression that does not require restricting the functional depen- 
dence of the bras and kets. 
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B. Second Look at the Nonstandard Expression 

In interpreting ( |72D we have concluded above that we must first regularize 
the integrand in order to obtain an integrable function. For the kinematic 
term J vr"^ Qab dh; dt — and even for the diffeomorphism constraint contribu- 
tion —J N°'Ha d^x dt — any natural regularization, such as one based on the 
expansion functions {hp{x)} and {/p^l^)}; ^ lattice formulation as dis- 
cussed in Sec. II, will be compatible regularizations. For these terms alone, 
the limit of the regularized functional integral as the regularization is re- 
moved will converge to the desired result. The implication of this fact is 
that these parts of the integrand are compatible with the initial (ultralocal) 
representation of the field operators; in fact, they are compatible with any 
diffeomorphism invariant realization. However, when it comes time to con- 
sider the Hamiltonian constraint, the behavior is quite different. While it is 
true that regularizing the Hamiltonian constraint will lead to a set of well- 
defined functional integrals, the limit of such regularized expressions will not 
converge to an acceptable result. There are two basic and important reasons 
for this unsatisfactory behavior, one of which (wrong field operator repre- 
sentation) we will deal with in this section, the other of which (perturbative 
nonrenormalizability of gravity) we will discuss in the next section. 

The first reason for the lack of a suitable convergence of the regularized 
nonstandard functional integral relates to the fact that the representation of 
the field operators needed to satisfy the Hamiltonian constraint is unitarily 
mequivalent to the ultralocal representation imposed in the initial stage of 
the quantization procedure. It is at the present stage of the analysis that we 
finally encounter the fact that our initial choice of field operator representa- 
tion is incompatible with making the Hamiltonian constraint operator T-i{x) 
into a densely defined local operator. Stated otherwise, using the ultralocal 
operator representation, the operator J N{x,t)T-C{x) d^x dt, for any nonzero 
smooth function N{x,t), has only the zero vector in its domain. This defect 
must be fixed before proceeding, and in so doing we will be explicitly led 
to a new representation of the field operators, one that is unitarily inequiv- 
alent to our starting (ultralocal) representation. In the process of effecting 
this change of representation, the scalar density b{x) will disappear from the 
scene entirely! 
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Pedagogical example 

It is pedagogically useful to outline an analogous story for a simpler and more 
familiar example. Consider general, locally self-adjoint field and momentum 
operators, (f){x) and n{x), x G M^, which satisfy the canonical commutation 
relations 

[^{x),n{y)]^iS{x-y). (78) 

Build a set of coherent states 

|7r,</.) = e*['^W-^('^)l|r7) , (79) 

where 0(7r) = J(f){x)7r{x) and 7r(0) = J7t{x)(f){x) (fy, with tt and (p real 
test functions, and \r]) is a normalized but otherwise unspecified fiducial 
vector. Note well that the choice of \ri) in effect determines the representation 
of the canonical field operators. We next present a portion of the story from 
Ref. 8. 

We initially choose I77) to correspond to an ultralocal representation such 
that 

(tt", (/.'V, (/.') = exp{ii/[(/."(x)7r'(a;) - tt" (x) <!)' (x)] d^x} 
X exp{-lJ{M{x)-'[n'\x) - n\x)f + M{x)[(P''{x) - (f)'ix)f} d'^x) . (80) 

Here M{x), < M{x) < 00, is an arbitrary (smooth) function of the ul- 
tralocal representation with the dimensions of M. The given ultralocal field 
operator representation is in fact unitarily inequivalent for each distinct func- 
tion M{x). [Note well that M{x) here plays the role of b{x) in the present 
paper.] 

We wish to apply this formulation to describe the relativistic free field of 
mass m for which the Hamiltonian operator is formally given by 

H^\j: {T^{xf + [V^{x)f + m^{xf} : dh . (81) 

If we build this operator out of the field and momentum operators in the 
ultralocal representation, then no matter what vector is used to define : :, 
V, will have only the zero vector in its domain. We need to change the 
field operator representation, which means we have to change the fiducial 
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vector from \r]) to |0; m), the true ground state of the proposed Hamiltonian 
operator 7i. 

Let us first regularize the formal Hamiltonian H. To that end, let 
denote a complete set of real, orthonormal functions on and define the 
sequence of kernels, for all N e {1, 2, 3, . . .}, given by 

N 

Kj^{x,y) = ^ Un{x)un{y) , (82) 

n=l 

which converges to S{x — y) as a distribution when N —>■ oo. Let 

M^) = jKj,{x,y)^{y)d\ , (83) 
t:n{x) = JKn(x, y)n(y) d\ , (84) 

and with these operators build the sequence of regularized Hamiltonian op- 
erators 

nN = U ■■ {7r7v(a;)2 + [VM^)]^ + m^M^)^} ■ d\ (85) 

for all N ^ where : : denotes normal order with respect to the ground state 
|0; m)jv of T-Ln- 

We would like to have a constructive way to identify the ground state of 
T-Ln- For this purpose consider the set 

= \ w — 1 —\ —\ : J <oo\ (86) 

for general sets {oj} (not all zero), {tTj}, and {0j}. (How these expressions 
may be generated is discussed in Ref. 8.) As grows, the general element 
in Sn becomes exponentially small, save for elements that correspond to 
vectors which well approximate the ground state |0;m)7v Suitable linear 
combinations can convert the original reproducing kernel (tt", 0"|7r', 0') to the 
reproducing kernel (tt", 0"; m|7r', 0'; m)7v which is based on a fiducial vector 
that has the form |0; m) n for the first N degrees of freedom and is unchanged 
for the remaining degrees of freedom. Finally, we may take the limit — > oo 
which then leads to 

(7r",0";m|7r',0';m) = exp{ltj[^"*{k)fr'{k)-n"*{k)^'{k)]d^ 
X exp{-\J{uj{k)-^\r{k) - 7r'{k)\^ + uj{k)\^"{k) - 4>'{k)\^} d^k) , (87) 
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where uj(k) = \/h^~+lrt? and n^k) = (27r)^^^'^ Je^^'^'^ n(x) d^x, etc. The pro- 
cedure sketched above is referred to as recentering the coherent states or 
equivalently as recentering the reproducing kernel. This form of reproducing 
kernel is no longer ultralocal and contains no trace of the scalar function 
M{x), whatever form it may have had. Moreover, and this is an important 
point, the new representation is fully compatible with the Hamiltonian Ti 
being a nonnegative, self-adjoint operator. Indeed, the expression for the 
propagator is given by 

(7r",0";m|e-^^^|7r',0';m) = L" L' exp[JCik)e-''^^'^'' ('(k) d\] , (88) 

where 

C{k) = [uj{ky/^4>{k) + tuj{k)-^/^n{k)]/V2 , (89) 

L^exp[-lJ\ak)\^d'k]. (90) 

The definition offered by (|88D is continuous in T, which is the principal 
guarantor that the expression 

1-L = \j : {7r(x)2 + [V0(x)]2 + m20(x)2} : dh , (91) 

where : : denotes normal ordering with respect to the ground state |0;m) of 
the operator ?i, is a self-adjoint operator as desired. 

Let us summarize the basic content of the present pedagogical exam- 
ple. Even though we started with a very general ultralocal representation, 
as characterized by the general function M{x), we have forced a complete 
change of representation to one compatible with the Hamiltonian operator 
for a relativistic free field of arbitrary mass m. In so doing all trace of the 
initial arbitrary function M{x) has disappeared, and in its place, effectively 
speaking, has appeared the pseudo-differential operator V— -|- m? having 
only its dimension (mass) in common with the original function M{x). The 
original ultralocal representation is completely gone! (Remark: A moments 
reflection should convince the reader that a comparable analysis can be made 
for either the interacting 02 03 model as well, both of which satisfy ([78|) , 
showing that the general argument is not limited just to free theories; see 
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Strong coupling gravity 



The discussion in the present paper has been predicated on the assumption 
that we are analyzing the gravitational field and therefore the classical Hamil- 
tonian is that given in (^Bj). However, it is pedagogically instructive if we 



briefiy comment on an approximate theory — based on the so-called "strong 
coupling approximation" — where the Hamiltonian constraint is re- 
placed by the expression 

HscA{x)^g{x)-'/^[K{x)n'a{^) - \K{x)t:1{x)] + 2Kg{xf/' , (92) 

in which the term proportional to '^'^'>R{x) has been dropped, and where we 
have also introduced the cosmological constant A (with dimensions L~^) and 
an associated auxiliary term in the Hamiltonian. The result is a model for 
which the new Hamiltonian constraint (^) is indeed compatible with some 
form of an ultralocal representation. The proper form of that ultralocal rep- 
resentation may be determined by a similar procedure, i.e., by studying an 
analogue of the set S'tv, and by recentering the reproducing kernel based on 
ensuring that the quantum version of Hsca{x) is a local self-adjoint opera- 
tor. In so doing, we note that it may happen that not all arbitrariness of the 
original scalar density b{x) is "squeezed out" by the recentering procedure 
described above. This situation occurs because a one-parameter arbitrari- 
ness generally remains for typical ultralocal models [0. (Any remaining 
arbitrariness in the ultralocal case is in contrast with that of the true Hamil- 
tonian constraint for gravity for which we expect no trace of the original 
function b{x) to remain.) These remarks conclude our discussion of strong 
coupling gravity. 



C. Third Look at the Nonstandard Expression 

The preceding discussion has been based on the assumption that some fiducial 
vector can be found compatible with the Hamiltonian operator constraint, or 
stated otherwise, that the Hamiltonian constraint can actually be realized as 
a local self-adjoint operator. This requirement is by no means obvious, and it 
is to this issue that we now turn our attention. The difficulty arises because 
the naive form of the Hamiltonian constraint operator almost surely needs 
some form of renormalization if it is going to be well defined. If perturbation 
theory is any guide, we not only expect that there will be renormalization 
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counterterms, but because gravity is perturbatively nonrenormalizable, one 
may expect an infinite number of distinct counterterms. On tlie otlier liand, 
as we next argue, it is possible tliat perturbation tlieory is not a very reliable 
guide in the case of perturbatively nonrenormalizable theories. 



Nonrenormalizable scalar fields 

Consider the case of perturbatively nonrenormalizable quartic, self-interact- 
ing scalar fields, i.e., the so-called 0^ theories, where the space-time dimen- 
sion n > 5. On the one hand, viewed perturbatively, such theories entail an 
infinite number of distinct counterterms. On the other hand, the continuum 
limit of a straightforward Euclidean lattice formulation leads to a quasifree 
theory — a genuinely noninteracting theory — whatever choice is made for the 
renormalized field strength, mass, and coupling constant p4[. In the author's 



view both of these results are unsatisfactory. Instead, it is possible that an 
intermediate behavior holds true, even though that cannot be proven yet. 
Let us illustrate an analogous but simpler situation where the conjectured 
intermediate behavior can be rigorously established. 

Consider an ultralocal quartic interacting scalar field, which, viewed clas- 
sically, is nothing but the relativistic 0^ model with all the spatial gradients 
in the usual free term dropped. As a mathematical model of quantum field 
theory, an ultralocal model is readily seen to be perturbatively nonrenor- 
malizable, while the continuum limit of a straightforward lattice formulation 
becomes quasifree, basically because of the vise grip of the Central Limit The- 
orem. Perturbative nonrenormalizability and lattice-limit triviality is similar 
to the behavior for relativistic 0^ models, but for the simpler ultralocal model 



for which an intermediate approach can be rigorously proven to hold [|I2 
Roughly speaking, a characterization of this intermediate behavior is the fol- 
lowing: From a functional integral standpoint, and for any positive value of 
the quartic coupling constant, the quartic interaction acts like a hard-core 
in history space projecting out certain contributions that would otherwise 
be allowed by the free theory alone. This phenomenon takes the form of a 
nonstandard, nonclassical counterterm in the Hamiltonian that does not van- 
ish as the coupling constant of the quartic interaction vanishes. Specifically, 
for the model in question, the additional counterterm is a counterterm for 
the kinetic energy and is formally proportional to ff/(f){x)'^, which in form 
is not unlike the centripetal potential that arises in spherical coordinates in 
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three-dimensional quantum mechanics. In summary, inclusion of a formal 
additional interaction proportional to in the Hamiltonian density is 

sufficient to result in a well defined and nontrivial (i.e., non-Gaussian) quan- 
tum theory for interacting ultralocal scalar models. In addition, it may be 
shown [T^ that the classical limit of such quantum theories reproduces the 



classical model with which one started. 

The foregoing brief summary holds rigorously for the ultralocal scalar 
fields, and it is conjectured that a suitable counterterm would lead to an ac- 
ceptable intermediate behavior for the relativistic models 0^, n > 5. What 
form should the counterterm take in the case of the relativistic 0^ mod- 
els? We can make a plausible suggestion guided by the following general 
principle that holds in the ultralocal case: The counterterm should be an 
ultralocal (because the kinetic energy is ultralocal) potential term arising 
from the kinetic energy. For the relativistic field that argument suggests the 
counterterm should again be proportional to h'^/(j){xY. It is also part of this 
general conjecture that the same counterterm is not limited to 0^ models, 
but should be effective for other nonrenormalizable interactions, e.g., such 
as 0^, n > 4, etc. The full argument available at present to support this 
conjecture appears in Chap. 11 of [Q. (It may even be possible to exam this 
proposal by means of suitable Monte Carlo studies, but so far this challenge 
has not been taken up.) 

Note well that the hard-core picture of nonrenormalizable interactions 
leads to such interactions behaving as discontinuous perturbations: Once 
turned on, such interactions cannot be completely turned off! Stated oth- 
erwise, as the nonlinear coupling constant is reduced toward zero, the the- 
ory passes continuously to a "pseudofree" theory — different than the "free" 
theory — which retains the effects of the hard core. The interacting theory is 
continuously connected to the pseudofree theory, and may even possess some 
form of perturbation theory about the pseudofree theory. Evidently, the pres- 
ence of the hard-core interaction makes any perturbation theory developed 
about the original unperturbed theory almost totally meaningless. 

Nonrenormalizable gravity 

Although the differences between gravity and nonrenormalizable scalar inter- 
action are significant in their details, there are certain similarities we wish to 



draw on. Most importantly, one can argue p5[ that the nonlinear contribu- 



37 



tions to gravity act as a hard-core interaction in a quantization scheme, and 
thus the general picture sketched above for nonrenormahzable scalar fields 
should apply to gravity as well. Assuming that the analogy holds further, 
there should be a nonstandard, nonclassical counterterm that incorporates 
the dominant, irremovable effects of the hard-core interaction. Accepting 
the principle that in such cases perturbation theory offers no clear hint as 
to what counterterms should be chosen, we appeal to the guide used in the 
scalar case. Thus, as our proposed counterterm, we look for an ultralocal 
potential arising from the kinetic energy that appears in the Hamiltonian 
constraint. In fact, the only ultralocal potential that has the right transfor- 
mation properties is proportional to fi? g^xY^"^. Thus we are led to conjecture 
that the "nonstandard counterterm" is none other than a term like the famil- 
iar cosmological constant contribution! Unlike the scalar field which required 
an unusual term proportional to the gravitational case has resulted 

in suggesting a term proportional to an "old friend", namely g^xY/"^. At 
first glance, it seems absurd that such a harmless looking term could act to 
"save" the nonrenormalizability of gravity. In its favor we simply note that 
the analogy with how other nonrenormahzable theories are "rescued" is too 
strong to dismiss the present proposal out of hand — and of course one must 
resist any temptation "to think perturbatively" . Any attempt to consider 
this possibility must wait until another occasion; we hope to return to this 
subject elsewhere. 

As one small aspect of this problem, let us briefly discuss how the factor 
arises in the gravitational case. Merely from a dimensional point of view, 

we note that (the first term of) the local kinetic energy operator has a formal 

structure given by 

-^^^^^i c \ . 9ac{x) gix)-^/"^ gbd{x) \ ) , (93) 

where we have restored the factor IGTrG/c^. Thus the anticipated countert- 
erm is proportional to {GY? j c')g{xf'l'^ . We next cast this term into the 
usual form for a contribution to the potential part of the Hamiltonian con- 
straint, namely, in a form proportional to {(? j G) Kg{xY^'^ ■ Hence, to recast 
our anticipated counterterm into this form, we need a factor proportional to 

^ - l%ian.. - (10-^^cm)^ . (94) 
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In the classical symbol for the Hamiltonian constraint operator, this factor is 
multiplied by an expectation value with dimensions originating from the 
density nature of the two momentum factors and leading to an overall factor 
with the dimensions that is proportional to as claimed. Let us call 
the resultant counterterm Ac gixY^"^- Since the sign of the DeWitt metric 
that governs the kinetic energy term is indefinite, it is not even possible 
to predict the sign of A^;. However, one thing appears certain. While the 
proposed counterterm Ac gixY^"^ is certainly not cosmological in origin, its 
influence may well be! 

The foregoing scenario has assumed the hard-core model of nonrenormal- 
izable interactions applies to the theory of gravity. However, that may well 
not be the case, and, instead, some other counterterm (s) may be required 
to cure the theory of gravity. Note well that the general structure of our 
approach to quantize gravity is largely insensitive to just what form of regu- 
larization and renormalization is required. In particular, the use of the affine 
field variables, the application of the projection operator method to impose 
constraints, and the development of the nonstandard phase-space functional 
integral representation for the reproducing kernel of the regularized physical 
Hilbert space all have validity quite independently of the form in which the 
Hamiltonian constraint is ultimately turned into a local sclf-adjoint operator. 
Although we have outlined one particular version in which the Hamiltonian 
constraint may possibly be made into a densely defined local operator, we are 
happy to keep an open mind about the procedure by which this ultimately 
may take place since many different ways in which this process can occur 
arc fully compatible with the general principles of our proposed quantization 
scheme for the gravitational field. 
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